Alzheimer's polygenic risk scores are associated with cognitive phenotypes in Down syndrome.

INTRODUCTION: This study aimed to investigate the influence of the overall Alzheimer's disease (AD) genetic architecture on Down syndrome (DS) status, cognitive measures, and cerebrospinal fluid (CSF) biomarkers. METHODS: AD polygenic risk scores (PRS) were tested for association with DS-related traits. RESULTS: The AD risk PRS was associated with disease status in several cohorts of sporadic late- and early-onset and familial late-onset AD, but not in familial early-onset AD or DS. On the other hand, lower DS Mental Status Examination memory scores were associated with higher PRS, independent of intellectual disability and APOE (PRS including APOE, PRSAPOE , p = 2.84 × 10-4 ; PRS excluding APOE, PRSnonAPOE , p = 1.60 × 10-2 ). PRSAPOE exhibited significant associations with Aß42, tTau, pTau, and Aß42/40 ratio in DS. DISCUSSION: These data indicate that the AD genetic architecture influences cognitive and CSF phenotypes in DS adults, supporting common pathways that influence memory decline in both traits. HIGHLIGHTS: Examination of the polygenic risk of AD in DS presented here is the first of its kind. AD PRS influences memory aspects in DS individuals, independently of APOE genotype. These results point to an overlap between the genes and pathways that leads to AD and those that influence dementia and memory decline in the DS population. APOE ε4 is linked to DS cognitive decline, expanding cognitive insights in adults.

Potential risk factors and genetic variants associated with dental caries incidence in Appalachia using genome-wide survival analysis.

OBJECTIVE: The aim of this study was to identify the potential risk factors and genetic variants associated with dental caries incidence using survival analysis. METHODS: The Center for Oral Health Research in Appalachia recruited and prospectively followed pregnant women and their children. A total of 909 children followed from birth for up to 7 years were included in this study. Annual intra-oral examinations were performed to assess dental caries experience including the approximate time to first caries incidence in the primary dentition. Cox proportional hazards models were used to assess the associations of time to first caries incidence with self-reported risk factors and 4.9 million genetic variants ascertained using a genome-wide genotyping array. RESULTS: A total of 196 of 909 children (21.56%) had their first primary tooth caries event during follow-up. Household income, home water source, and mother's educational attainment were significantly associated with time to first caries incidence in the stepwise Cox model. The heritability (i.e., proportion of variance explained by genetics) of time to first caries was 0.54. Though no specific genetic variants were associated at the genome-wide significance level (P < 5E-8), we identified 14 loci at the suggestive significance level (5E-8 < P < 1E-5), some of which were located within or near genes with plausible biological functions in dental caries. CONCLUSION: Our findings indicate that household income, home water source, and mother's educational attainment are independent risk factors for dental caries incidence. We nominate several suggestive loci for further investigation.

Genome-wide analysis identifies novel loci influencing plasma apolipoprotein E concentration and Alzheimer's disease risk.

The APOE 2/3/4 polymorphism is the greatest genetic risk factor for Alzheimer's disease (AD). This polymorphism is also associated with variation in plasma ApoE level; while APOE*4 lowers, APOE*2 increases ApoE level. Lower plasma ApoE level has also been suggested to be a risk factor for incident dementia. To our knowledge, no large genome-wide association study (GWAS) has been reported on plasma ApoE level. This study aimed to identify new genetic variants affecting plasma ApoE level as well as to test if baseline ApoE level is associated with cognitive function and incident dementia in a longitudinally followed cohort of the Ginkgo Evaluation of Memory (GEM) study. Baseline plasma ApoE concentration was measured in 3031 participants (95.4% European Americans (EAs)). GWAS analysis was performed on 2580 self-identified EAs where both genotype and plasma ApoE data were available. Lower ApoE concentration was associated with worse cognitive function, but not with incident dementia. As expected, the risk for AD increased from E2/2 through to E4/4 genotypes (P for trend = 4.8E-75). In addition to confirming the expected and opposite associations of APOE*2 (P = 4.73E-79) and APOE*4 (P = 8.73E-12) with ApoE level, GWAS analysis revealed nine additional independent signals in the APOE region, and together they explained about 22% of the variance in plasma ApoE level. We also identified seven new loci on chromosomes 1, 4, 5, 7, 11, 12 and 20 (P range = 5.49E-08 to 5.36E-10) that explained about 9% of the variance in ApoE level. Plasma ApoE level-associated independent variants, especially in the APOE region, were also associated with AD risk and amyloid deposition in the brain, indicating that genetically determined ApoE level variation may be a risk factor for developing AD. These results improve our understanding of the genetic determinants of plasma ApoE level and their potential value in affecting AD risk.

DNA methylation and 28-year cardiovascular disease risk in type 1 diabetes: the Epidemiology of Diabetes Complications (EDC) cohort study.

BACKGROUND: The potential for DNA methylation (DNAm) as an early marker for cardiovascular disease (CVD) and how such an association might differ by glycemic exposure has not been examined in type 1 diabetes, a population at increased CVD risk. We thus performed a prospective epigenome-wide association study of blood leukocyte DNAm (EPIC array) and time to CVD incidence over 28 years in a childhood-onset (< 17 years) type 1 diabetes cohort, the Pittsburgh Epidemiology of Diabetes Complications (EDC) study (n = 368 with DNA and no CVD at baseline), both overall and separately by glycemic exposure, as measured by HbA1c at baseline (split at the median: < 8.9% and ≥ 8.9%). We also assessed whether DNAm-CVD associations were independent of established cardiometabolic risk factors, including body mass index, estimated glucose disposal rate, cholesterol, triglycerides, blood pressure, pulse rate, albumin excretion rate, and estimated glomerular filtration rate. RESULTS: CVD (first instance of CVD death, myocardial infarction, coronary revascularization, ischemic ECG, angina, or stroke) developed in 172 participants (46.7%) over 28 years. Overall, in Cox regression models for time to CVD, none of the 683,597 CpGs examined reached significance at a false discovery rate (FDR) ≤ 0.05. In participants with HbA1c < 8.9% (n = 180), again none reached FDR ≤ 0.05, but three were associated at the a priori nominal significance level FDR ≤ 0.10: cg07147033 in MIB2, cg12324048 (intergenic, chromosome 3), and cg15883830 (intergenic, chromosome 1). In participants with HbA1c ≥ 8.9% (n = 188), two CpGs in loci involved in calcium channel activity were significantly associated with CVD (FDR ≤ 0.05): cg21823999 in GPM6A and cg23621817 in CHRNA9; four additional CpGs were nominally associated (FDR ≤ 0.10). In participants with HbA1c ≥ 8.9%, DNAm-CVD associations were only modestly attenuated after cardiometabolic risk factor adjustment, while attenuation was greater in those with HbA1c < 8.9%. No pathways were enriched in those with HbA1c < 8.9%, while pathways for calcium channel activity and integral component of synaptic membrane were significantly enriched in those with HbA1c ≥ 8.9%. CONCLUSIONS: These results provide novel evidence that DNAm at loci involved in calcium channel activity and development may contribute to long-term CVD risk beyond known risk factors in type 1 diabetes, particularly in individuals with greater glycemic exposure, warranting further study.

Rare variants found in clinical gene panels illuminate the genetic and allelic architecture of orofacial clefting.

PURPOSE: Orofacial clefts (OFCs) are common birth defects including cleft lip, cleft lip and palate, and cleft palate. OFCs have heterogeneous etiologies, complicating clinical diagnostics because it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs; therefore, we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls. METHODS: We evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria. RESULTS: 9.04% of cases and 1.02% of controls had "likely pathogenic" variants (P < .0001), which was almost exclusively driven by heterozygous variants in autosomal genes. Cleft palate (17.6%) and cleft lip and palate (9.09%) cases had the highest yield, whereas cleft lip cases had a 2.80% yield. Out of 39 genes with likely pathogenic variants, 9 genes, including CTNND1 and IRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were "variants of uncertain significance", occurring more frequently in cases (P = .004), but no individual gene showed a significant excess of variants of uncertain significance. CONCLUSION: These results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

Meta-analysis of age-related cognitive decline reveals a novel locus for the attention domain and implicates a COVID-19-related gene for global cognitive function.

INTRODUCTION: Cognitive abilities have substantial heritability throughout life, as shown by twin- and population-based studies. However, there is limited understanding of the genetic factors related to cognitive decline in aging across neurocognitive domains. METHODS: We conducted a meta-analysis on 3045 individuals aged ≥65, derived from three population-based cohorts, to identify genetic variants associated with the decline of five neurocognitive domains (attention, memory, executive function, language, visuospatial function) and global cognitive decline. We also conducted gene-based and functional bioinformatics analyses. RESULTS: Apolipoprotein E (APOE)4 was significantly associated with decline of memory (p = 5.58E-09) and global cognitive function (p = 1.84E-08). We identified a novel association with attention decline on chromosome 9, rs6559700 (p = 2.69E-08), near RASEF. Gene-based analysis also identified a novel gene, TMPRSS11D, involved in the activation of SARS-CoV-2, to be associated with the decline in global cognitive function (p = 4.28E-07). DISCUSSION: Domain-specific genetic studies can aid in the identification of novel genes and pathways associated with decline across neurocognitive domains. HIGHLIGHTS: rs6559700 was associated with decline of attention. APOE4 was associated with decline of memory and global cognitive decline. TMPRSS11D, a gene involved in the activation of SARS-CoV-2, was implicated in global cognitive decline. Cognitive domain abilities had both unique and shared molecular pathways across the domains.

Genomic landscape of Down syndrome-associated acute lymphoblastic leukemia.

Trisomy 21, the genetic cause of Down syndrome (DS), is the most common congenital chromosomal anomaly. It is associated with a 20-fold increased risk of acute lymphoblastic leukemia (ALL) during childhood and results in distinctive leukemia biology. To comprehensively define the genomic landscape of DS-ALL, we performed whole-genome sequencing and whole-transcriptome sequencing (RNA-Seq) on 295 cases. Our integrated genomic analyses identified 15 molecular subtypes of DS-ALL, with marked enrichment of CRLF2-r, IGH::IGF2BP1, and C/EBP altered (C/EBPalt) subtypes compared with 2257 non-DS-ALL cases. We observed abnormal activation of the CEBPD, CEBPA, and CEBPE genes in 10.5% of DS-ALL cases via a variety of genomic mechanisms, including chromosomal rearrangements and noncoding mutations leading to enhancer hijacking. A total of 42.3% of C/EBP-activated DS-ALL also have concomitant FLT3 point mutations or insertions/deletions, compared with 4.1% in other subtypes. CEBPD overexpression enhanced the differentiation of mouse hematopoietic progenitor cells into pro-B cells in vitro, particularly in a DS genetic background. Notably, recombination-activating gene-mediated somatic genomic abnormalities were common in DS-ALL, accounting for a median of 27.5% of structural alterations, compared with 7.7% in non-DS-ALL. Unsupervised hierarchical clustering analyses of CRLF2-rearranged DS-ALL identified substantial heterogeneity within this group, with the BCR::ABL1-like subset linked to an inferior event-free survival, even after adjusting for known clinical risk factors. These results provide important insights into the biology of DS-ALL and point to opportunities for targeted therapy and treatment individualization.

Investigation of the independent role of a rare APOE variant (L28P; APOE*4Pittsburgh) in late-onset Alzheimer disease.

A rare missense APOE variant (L28P; APOE*4Pittsburgh), which is present only in populations with European ancestry, has been reported to be a risk factor for late-onset Alzheimer's disease (LOAD). However, due to the complete linkage disequilibrium of L28P with APOE*4 (C112R), its independent genetic association is uncertain. The original association study implicating L28P with LOAD risk was carried out in a relatively small sample size. In the current study, we have re-evaluated this association in a large case-control sample of 15,762 White U.S. subjects and investigated its independent effect in APOE 3/4 subjects, as L28P has been observed only in the heterozygous state of APOE*4 carriers and 3/4 is the most common genotype containing the APOE*4 allele. The heterozygous carrier frequency of L28P, all with APOE*4, was about 3-fold higher in AD cases than in cognitively intact controls (0.845% vs. 0.277%). The age- and sex-adjusted meta-analysis odds ratio (OR) was 2.87 (95% CI: 1.34 - 6.13; = 0.0066). Among APOE 3/4 subjects, age- and sex-adjusted meta-analysis OR was 1.53 (95% CI: 0.70 - 3.36; p = 0.28), indicating its effect was independent of APOE*4. The lack of statistical significance appears mainly due to the low power of 4138 subjects with the 3/4 genotype (12% power at α= 0.05) compared to the required sample of 139,088 subjects with the 3/4 genotype to detect an OR of 1.5 at α= 0.05 and 80% power. Our data suggesting that L28P has an independent genetic effect on AD risk is reinforced by earlier experimental findings showing that this mutation leads to significant structural and conformational changes in the ApoE4 molecule and can induce functional defects associated with neuronal Aß42 accumulation and oxidative stress. Additional functional studies in cell-based systems and animal models will help to delineate its functional significance in the etiology of AD.

Association of Fc Gamma Receptor 3B Gene Copy Number Variation with Rheumatoid Arthritis Susceptibility.

Structural variations such as copy number variants (CNVs) have been associated with multiple autoimmune diseases. In this study, we explored the association of the Fc gamma receptor 3B gene (FCGR3B) copy number variation (CNV) with rheumatoid arthritis (RA) susceptibility and related serological traits in the Pakistani population. We also performed a meta-analysis of four published FCGR3B CNV studies along with the current study. A total of 927 subjects (597 RA cases, 330 healthy controls) were recruited from three rheumatology centers in Pakistan. Anti-cyclic citrullinated peptide (anti-CCP) antibodies and rheumatoid factor (RF) were measured in RA patients. FCGR3B copy number was assayed using the TaqMan® CN assay (Hs04211858_cn, Applied Biosystems, Foster City, CA, USA) and the copy number was estimated by using CopyCaller® software (version 2.1; Applied Biosystems, USA). Logistic regression was applied to calculate the odds ratio (OR) of RA risk associated with FCGR3B CNV using sex and age as covariates in R. Meta-analysis on four previously published studies and the current study was performed using the random-effect model. We observed a significant association between FCGR3B copy number < 2 and RA susceptibility (OR = 1.53; 95% CI: 1.05 to 2.22; p = 0.0259) and anti-CCP seropositivity (OR 2.56; 95% CI: 1.34 to 4.89; p = 0.0045). A non-significant association of FCGR3B copy number < 2 was also observed between increased rheumatoid factor (RF) seropositivity (OR = 1.74; 95% CI:0.93 to 3.26; p = 0.0816). Meta-analysis on 13,915 subjects (7005 RA cases and 6907 controls) also showed significant association of copy number < 2 with the increased risk of RA (OR = 1.30; 95% CI: 1.07 to 1.56; p = 0.00671). FCGR3B copy number < 2 is associated with increased RA risk and anti-CCP seropositivity.

Genome-Wide Association Study of Incident Dementia in a Community-Based Sample of Older Subjects.

BACKGROUND: Alzheimer's disease (AD) is a complex disease influenced by the environment and genetics; however, much of the genetic component remains unaccounted for. OBJECTIVE: The purpose of this work was to use genome-wide association analyses to detect genetic associations with incident AD in a sample of older adults aged 75 and above. METHODS: We performed a genome-wide association study (GWAS) on genome-wide genotyped and imputed data (14,072,053 variants) on the Gingko Evaluation of Memory (GEM) study sample consisting of 424 incident dementia (mean age = 84.46±3.91) and 2,206 non-demented (mean age = 84.55±3.23) subjects. RESULTS: The established association of APOE*4 carriers with AD was confirmed in this community-based sample of older subjects (odds ratio (OR) = 2.22; p = 9.36E-14) and was stronger in females (OR = 2.72; p = 1.74E-10) than in males (OR = 1.88; p = 2.43E-05). We observed a novel genome-wide significant (GWS) locus on chromosome 12 near ncRNA LOC105369711/rs148377161 (OR = 3.31; p = 1.66E-08). In addition, sex-stratified analyses identified two novel associations in males: one near ncRNA LOC729987/rs140076909 on chromosome 1 (OR = 4.51; p = 3.72E-08) and the other approaching GWS near ncRNA LOC105375138/rs117803234 on chromosome 7 (OR = 3.76; p = 6.93E-08). CONCLUSION: The use of community-based samples of older individuals and incident dementia as a phenotype may be a helpful approach for the identification of novel genes for AD, which may not be detected in standard case-control studies. Replication of these signals and further studies of these regions and genes will help to provide a clearer picture for their role in AD.

Genome-wide association study of multiethnic nonsyndromic orofacial cleft families identifies novel loci specific to family and phenotypic subtypes.

Nonsyndromic orofacial clefts (OFCs) are among the most common craniofacial birth defects worldwide, and known to exhibit phenotypic and genetic heterogeneity. Cleft lip plus cleft palate (CLP) and cleft lip only (CL) are commonly combined together as one phenotype (CL/P), separately from cleft palate alone. In comparison, our study analyzes CL and CLP separately. A sample of 2218 CL and CLP cases, 4537 unaffected relatives of cases, and 2673 pure controls with no family history of OFC were selected from the Pittsburgh Orofacial Cleft (Pitt-OFC) multiethnic study.genome-wide association studies were run for seven specific phenotypes created based on the cleft type(s) observed within these families, as well as the combined CL/P phenotype. Five novel genome-wide significant associations, 3q29 (rs62284390), 5p13.2 (rs609659), 7q22.1 (rs6465810), 19p13.3 (rs628271), and 20q13.33 (rs2427238), and nine associations (p ≤ 1.0E-05) within previously confirmed OFC loci-PAX7, IRF6, FAM49A, DCAF4L2, 8q24.21, ARID3B, NTN1, TANC2 and the WNT9B:WNT3 gene cluster-were observed. We also found that single nucleotide polymorphisms within a subset of the associated loci, both previously known and novel, differ substantially in terms of their effects across cleft- or family-specific phenotypes, indicating not only etiologic differences between CL and CLP, but also genetic heterogeneity within each of the two OFC subtypes.

Genome-wide Interaction Study Implicates VGLL2 and Alcohol Exposure and PRL and Smoking in Orofacial Cleft Risk.

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect, affecting approximately 1 in 700 births. NSCL/P has complex etiology including several known genes and environmental factors; however, known genetic risk variants only account for a small fraction of the heritability of NSCL/P. It is commonly suggested that gene-by-environment (G×E) interactions may help explain some of the "missing" heritability of NSCL/P. We conducted a genome-wide G×E interaction study in cases and controls of European ancestry with three common maternal exposures during pregnancy: alcohol, smoking, and vitamin use using a two-stage design. After selecting 127 loci with suggestive 2df tests for gene and G x E effects, 40 loci showed significant G x E effects after correcting for multiple tests. Notable interactions included SNPs of 6q22 near VGLL2 with alcohol and 6p22.3 near PRL with smoking. These interactions could provide new insights into the etiology of CL/P and new opportunities to modify risk through behavioral changes.

Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

Association of Early Childhood Caries with Bitter Taste Receptors: A Meta-Analysis of Genome-Wide Association Studies and Transcriptome-Wide Association Study.

Although genetics affects early childhood caries (ECC) risk, few studies have focused on finding its specific genetic determinants. Here, we performed genome-wide association studies (GWAS) in five cohorts of children (aged up to 5 years, total N = 2974, cohorts: Center for Oral Health Research in Appalachia cohorts one and two [COHRA1, COHRA2], Iowa Fluoride Study, Iowa Head Start, Avon Longitudinal Study of Parents and Children [ALSPAC]) aiming to identify genes with potential roles in ECC biology. We meta-analyzed the GWASs testing ~3.9 million genetic variants and found suggestive evidence for association at genetic regions previously associated with caries in primary and permanent dentition, including the ß-defensin anti-microbial proteins. We then integrated the meta-analysis results with gene expression data in a transcriptome-wide association study (TWAS). This approach identified four genes whose genetically predicted expression was associated with ECC (p-values < 3.09 × 10−6; CDH17, TAS2R43, SMIM10L1, TAS2R14). Some of the strongest associations were with genes encoding members of the bitter taste receptor family (TAS2R); other members of this family have previously been associated with caries. Of note, we identified the receptor encoded by TAS2R14, which stimulates innate immunity and anti-microbial defense in response to molecules released by the cariogenic bacteria, Streptococcus mutans and Staphylococcus aureus. These findings provide insight into ECC genetic architecture, underscore the importance of host-microbial interaction in caries risk, and identify novel risk genes.

PRICKLE1 × FOCAD Interaction Revealed by Genome-Wide vQTL Analysis of Human Facial Traits.

The human face is a highly complex and variable structure resulting from the intricate coordination of numerous genetic and non-genetic factors. Hundreds of genomic loci impacting quantitative facial features have been identified. While these associations have been shown to influence morphology by altering the mean size and shape of facial measures, their effect on trait variance remains unclear. We conducted a genome-wide association analysis for the variance of 20 quantitative facial measurements in 2,447 European individuals and identified several suggestive variance quantitative trait loci (vQTLs). These vQTLs guided us to conduct an efficient search for gene-by-gene (G × G) interactions, which uncovered an interaction between PRICKLE1 and FOCAD affecting cranial base width. We replicated this G × G interaction signal at the locus level in an additional 5,128 Korean individuals. We used the hypomorphic Prickle1 Beetlejuice (Prickle1 Bj ) mouse line to directly test the function of Prickle1 on the cranial base and observed wider cranial bases in Prickle1 Bj/Bj . Importantly, we observed that the Prickle1 and Focadhesin proteins co-localize in murine cranial base chondrocytes, and this co-localization is abnormal in the Prickle1 Bj/Bj mutants. Taken together, our findings uncovered a novel G × G interaction effect in humans with strong support from both epidemiological and molecular studies. These results highlight the potential of studying measures of phenotypic variability in gene mapping studies of facial morphology.

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations.

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10-8) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10-10). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.

FAT4 identified as a potential modifier of orofacial cleft laterality.

Orofacial clefts (OFCs) are common (1 in 700 births) congenital malformations that include a cleft lip (CL) and cleft lip and palate (CLP). These OFC subtypes are also heterogeneous themselves, with the CL occurring on the left, right, or both sides of the upper lip. Unilateral CL and CLP have a 2:1 bias towards left-sided clefts, suggesting a nonrandom process. Here, we performed a study of left- and right-sided clefts within the CL and CLP subtypes to better understand the genetic factors controlling cleft laterality. We conducted genome-wide modifier analyses by comparing cases that had right unilateral CL (RCL; N = 130), left unilateral CL (LCL; N = 216), right unilateral CLP (RCLP; N = 416), or left unilateral CLP (LCLP; N = 638), and identified a candidate region on 4q28, 400 kb downstream from FAT4, that approached genome-wide significance for LCL versus RCL (p = 8.4 × 10-8 ). Consistent with its potential involvement as a genetic modifier of CL, we found that Fat4 exhibits a specific domain of expression in the mesenchyme of the medial nasal processes that form the median upper lip. Overall, these results suggest that the epidemiological similarities in left- to right-sided clefts in CL and CLP are not reflected in the genetic association results.

Detecting Gene-Environment Interaction for Maternal Exposures Using Case-Parent Trios Ascertained Through a Case With Non-Syndromic Orofacial Cleft.

Two large studies of case-parent trios ascertained through a proband with a non-syndromic orofacial cleft (OFC, which includes cleft lip and palate, cleft lip alone, or cleft palate alone) were used to test for possible gene-environment (G × E) interaction between genome-wide markers (both observed and imputed) and self-reported maternal exposure to smoking, alcohol consumption, and multivitamin supplementation during pregnancy. The parent studies were as follows: GENEVA, which included 1,939 case-parent trios recruited largely through treatment centers in Europe, the United States, and Asia, and 1,443 case-parent trios from the Pittsburgh Orofacial Cleft Study (POFC) also ascertained through a proband with an OFC including three major racial/ethnic groups (European, Asian, and Latin American). Exposure rates to these environmental risk factors (maternal smoking, alcohol consumption, and multivitamin supplementation) varied across studies and among racial/ethnic groups, creating substantial differences in power to detect G × E interaction, but the trio design should minimize spurious results due to population stratification. The GENEVA and POFC studies were analyzed separately, and a meta-analysis was conducted across both studies to test for G × E interaction using the 2 df test of gene and G × E interaction and the 1 df test for G × E interaction alone. The 2 df test confirmed effects for several recognized risk genes, suggesting modest G × E effects. This analysis did reveal suggestive evidence for G × Vitamin interaction for CASP9 on 1p36 located about 3 Mb from PAX7, a recognized risk gene. Several regions gave suggestive evidence of G × E interaction in the 1 df test. For example, for G × Smoking interaction, the 1 df test suggested markers in MUSK on 9q31.3 from meta-analysis. Markers near SLCO3A1 also showed suggestive evidence in the 1 df test for G × Alcohol interaction, and rs41117 near RETREG1 (a.k.a. FAM134B) also gave suggestive significance in the meta-analysis of the 1 df test for G × Vitamin interaction. While it remains quite difficult to obtain definitive evidence for G × E interaction in genome-wide studies, perhaps due to small effect sizes of individual genes combined with low exposure rates, this analysis of two large case-parent trio studies argues for considering possible G × E interaction in any comprehensive study of complex and heterogeneous disorders such as OFC.

Genome-Wide Association Study of Non-syndromic Orofacial Clefts in a Multiethnic Sample of Families and Controls Identifies Novel Regions.

Orofacial clefts (OFCs) are among the most prevalent craniofacial birth defects worldwide and create a significant public health burden. The majority of OFCs are non-syndromic and vary in prevalence by ethnicity. Africans have the lowest prevalence of OFCs (~ 1/2,500), Asians have the highest prevalence (~1/500), Europeans and Latin Americans lie somewhere in the middle (~1/800 and 1/900, respectively). Thus, ethnicity appears to be a major determinant of the risk of developing OFC. The Pittsburgh Orofacial Clefts Multiethnic study was designed to explore this ethnic variance, comprising a large number of families and individuals (~12,000 individuals) from multiple populations worldwide: US and Europe, Asians, mixed Native American/Caucasians, and Africans. In this current study, we analyzed 2,915 OFC cases, 6,044 unaffected individuals related to the OFC cases, and 2,685 controls with no personal or family history of OFC. Participants were grouped by their ancestry into African, Asian, European, and Central and South American subsets, and genome-wide association run on the combined sample as well as the four ancestry-based groups. We observed 22 associations to cleft lip with or without cleft palate at 18 distinct loci with p-values < 1e-06, including 10 with genome-wide significance (<5e-08), in the combined sample and within ancestry groups. Three loci - 2p12 (rs62164740, p = 6.27e-07), 10q22.2 (rs150952246, p = 3.14e-07), and 10q24.32 (rs118107597, p = 8.21e-07) are novel. Nine were in or near known OFC loci - PAX7, IRF6, FAM49A, DCAF4L2, 8q24.21, NTN1, WNT3-WNT9B, TANC2, and RHPN2. The majority of the associations were observed only in the combined sample, European, and Central and South American groups. We investigated whether the observed differences in association strength were (a) purely due to sample sizes, (b) due to systematic allele frequency difference at the population level, or (c) due to the fact certain OFC-causing variants confer different amounts of risk depending on ancestral origin, by comparing effect sizes to observed allele frequencies of the effect allele in our ancestry-based groups. While some of the associations differ due to systematic differences in allele frequencies between groups, others show variation in effect size despite similar frequencies across ancestry groups.

The PAX1 locus at 20p11 is a potential genetic modifier for bilateral cleft lip.

Nonsyndromic orofacial clefts (OFCs) are a common birth defect and are phenotypically heterogenous in the structure affected by the cleft - cleft lip (CL) and cleft lip and palate (CLP) - as well as other features, such as the severity of the cleft. Here, we focus on bilateral and unilateral clefts as one dimension of OFC severity, because the genetic architecture of these subtypes is not well understood. We tested for subtype-specific genetic associations in 44 bilateral CL (BCL) cases, 434 unilateral CL (UCL) cases, 530 bilateral CLP cases (BCLP), 1123 unilateral CLP (UCLP) cases, and unrelated controls (N = 1626), using a mixed-model approach. While no novel loci were found, the genetic architecture of UCL was distinct compared to BCL, with 44.03% of suggestive loci having different effects between the two subtypes. To further understand the subtype-specific genetic risk factors, we performed a genome-wide scan for modifiers and found a significant modifier locus on 20p11 (p=7.53×10-9), 300kb downstream of PAX1, that associated with higher odds of BCL vs. UCL, and replicated in an independent cohort (p=0.0018) with no effect in BCLP (p>0.05). We further found that this locus was associated with normal human nasal shape. Taken together, these results suggest bilateral and unilateral clefts may have different genetic architectures. Moreover, our results suggest BCL, the rarest form of OFC, may be genetically distinct from the other OFC subtypes. This expands our understanding of modifiers for OFC subtypes and further elucidates the genetic mechanisms behind the phenotypic heterogeneity in OFCs.